Electric Field Effects on the Optical Vibrations in AB-Stacked Bilayer Graphene 
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Using first-principles methods, we show that an applied perpendicular electric field E breaks the 
inversion symmetry of AB-stacked bilayer graphene (BLG), thereby slightly mixing the two in-plane 
high-energy optical vibrations (E g and E u modes). The mixed amplitudes increase parabolically 
with respect to the field strength when E<2.0 V/nm, and then exhibit linear dependence when 
E>2.Q V/nm. In contrast, the mixing effect on the out-of-plane vibrations (Ai g and A^ u modes) 
is found to be much stronger, with the mixed amplitudes nearly an order of magnitude larger than 
those for the in-plane modes. For the two in-plane modes, we then calculate their phonon linewidths 
and frequency shifts as a function of the electric field as well as the Fermi level. Our results reveal 
delicate interplay between electrons and phonons in BLG, tunable by the applied fields and charge 
carrier densities. 
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AB-Stacked Bilayer graphene (BLG) is a unique plat- 
form that both the charge carrier densities (i.e., the Fermi 
level Ep) and the band structure can be tuned through 
applied dual-gate electric fields HHH. With an applied 
external electric field (EEF) of several V/nm, an energy 
gap opening of up to 250 meV has been reported in gated 
BLG [2]. Recently, the development of the electrolytic 
gate allows the charge carrier densities in BLG doping as 
high as \n\ « 2.4 x 10 13 cm -2 The tunability on the 
band structure and the Fermi level reveals rather intrigu- 
ing properties in this system, including the unusual high 
field transport [5[ , the renormalization of the phonon en- 
er gy [SM sl, and the Fano resonance in infrared spectra 
[si flo| - Interesting device application such as hot elec- 
tron bolometer has also been reported. 

Understanding the response of material's physical 
properties to external perturbations is crucial for its prac- 
tical applications. In particular, EEF affects both elec- 
tronic structure and lattice dynamics, and plays a key 
role on the performance of electronic devices. Although 
much attention has been paid to the field effects on the 
electron ic p roperties of BLG [H, 0, 0, [lH, 13] , very few 
studies [14j, LL5j have been conducted on the response of 
the lattice dynamics to the EEF. Using analytical meth- 
ods, Ando and Koshino [Hj studied the effects of an un- 
even charge doping and external electric fields on the self 
energy of the in-plane modes in BLG. On the other hand, 
it has been proposed more than forty years ago that EEF 
could be applied to break the symmetry of specific silent 
mode in SrTiOa so that it becomes Raman active jlfjl - 
flit . The electric field induced Raman scattering allows 
study of infrared-active and silent modes in crystals of 
high symmetry, thus probing the electronic properties of 
the system. Recent development of the tip enhanced Ra- 
man spectroscopy [19] provides controlled local probe of 
the response under a nonuniform electric field. However, 
it is unclear whether such an external perturbation will 



alter the phonon mode eigenvector itself or not, and to 
what extent the phonon wave function will be changed, 
if any. 

Raman and infrared (IR) spectroscopy arepowerful 
tools to probe the lattice dynamics of BLG 0, H, [H, [HI, 
2ol-22|. A recent IR measurement of gated BLG identi- 
fied large phonon linewidth (~30 cm -1 ) for the mode of 
frequency uo = 1587 cm -1 , which has been ascribed to 
the in-plane E u mode [9] , while another IR work [lOj as- 
cribed the observed mode to the E g mode because of bro- 
ken symmetry. A possible optical mode mixing induced 
by EEF is also proposed to account for the splitting of 
the G-band in the Raman measurement in BLG [2(1 l2lj. 
Clearly, a correct explanation of the experimental ob- 
servations requires a detailed knowledge of the response 
of these phonon modes to the EEF from first-principles. 
Such knowledge would also provide a necessary insight 
into the phonon, electron-phonon and their effects on the 
transport performance of device made of BLG 0, H3- 

In this work, we present a first-principles study of the 
electric field effects on the optical vibrations in BLG. 
The optical phonon modes with q = are calculated us- 
ing the density-functional perturbation theory (DFPT) 

251 as implemented in the Quantum ESPRESSO code 
[26[ at the local density approximation (LDA) level. 
Norm-conserving pseudopotential [27j for carbon has 
been adopted to describe the core- valence interactions. 
The wave functions of the valence electrons are expanded 
in plane waves with a kinetic energy cutoff of 70 Ry. A 
vacuum region of 20 A has been introduced to eliminate 
the artificial interaction between neighboring supercells 
along the z direction. The relaxed C-C bond length is 
1.42 A and the interlayer distance is 3.32 A for the 
BLG. Details have been presented in our previous work 

Hi Hj]- To mimic the electric field, a sawtooth potential 
is applied to the bilayer graphene sheets perpendicularly. 
More realistic configurations, such as substrate-induced 
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FIG. 1: (Color online) Schematic plots of the in-plane (a) high-branch mode with ujh — 1592 cm -1 and (b) low-branch mode 
with cjl=1587 cm -1 under zero and finite electric fields, respectively. The calculated modes £l(E) and £h(E) under field can 
be decomposed using £z,(0) and £h(0) as a complete basis set. The amplitude of (c) |cli,| 2 (\chh\ 2 ) and (d) |chx| 2 (\clh\ 2 ) 
changes with respect to the field. The dashed lines are parabolic fitting of the data as E < 2 V/nm. 




Electric Field (V/nm) 



FIG. 2: (Color online) Schematic plots of the out-of-plane (a) 
high-branch Am mode with uo — 894 cm -1 and (b) low-branch 
A\ g mode with cj — 891 cm -1 under zero electric field. The 
calculated modes £l(E) and ^h(E) under field can be decom- 
posed using £l(0) and £h(0) as a complete basis set. The am- 
plitude of (c) \c L l\ 2 (\chh\ 2 ) and (d) \c H l\ 2 (\c L h\ 2 ) changes 
with respect to the field. The dashed lines are parabolic fit- 
ting of the data as E < 2 V/nm. 



electric field [30], should not change the physics signifi- 
cantly because of similar band gap opening. Dipole cor- 
rections have been considered throughout all calculations 
below. 

Without EEF, BLG possesses a point group symmetry 
of Dsd- The two in-plane optical modes at q = are low- 



branch Eg (lj = 1587 cm -1 ) and high-branch E u (lj = 
1592 cm -1 ) modes, as schematically shown in Figs. 1(b) 
and 1(a), respectively. Depicted in Figs. 2(b) and 2(a) are 
the two out-of-plane optical modes A\ g (lj = 891 cm -1 ) 
and A^ u (lj = 894 cm -1 ) modes, respectively. We have 
calculated the phonon frequencies under various finite 
fields using DFPT. The variations of uj are within 2 cm -1 
for electric field up to 5 V/nm, indicating negligible field 
effects on the phonon frequency shift within the adiabatic 
approximation, in agreement with previous work [32| . 

We then consider the phonon mode mixing induced by 
EEF. The applied EEF perpendicular to BLG will break 
the inversion symmetry of BLG, and the point group is 
reduced to Cs v . The two in-plane modes now belong to 
the E representation of C% Vl and both of them become 
Raman and IR active. Since E is a two-dimensional rep- 
resentation, the two orthonormal eigenvectors of these 
modes would form a complete set for the mixed in-plane 
modes. Because the frequencies of the E g and E u modes 
are close to each other, a finite electric field will be very 
likely to alter the phonon wavefunctions. Hereafter the 
(mixed) E g and E u modes are denoted as low-branch 
(L) and high-branch (H) modes, respectively, as shown 
in Fig. 1. 

Quantitatively, the mixed amplitudes can be obtained 
by decomposition: £h(E) = c H l£,l(0) + c H h£,h(0), 
and £l(E) = c LL ^ L (^) + c L h£,h(0), with c a p = 
&(^)&(0) (a,j8 = L,H). The eigenvectors &(0) = 
[0.5, -0.5, 0.5, -0.5] T , and f# (0) = [0.5, -0.5, -0.5, 0.5] T 
correspond to the E g and E u modes under zero elec- 
tric field, respectively. The expansion coefficients sat- 
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isfy \c H l\ 2 + \c L l\ 2 = \c L h\ 2 + \c H h\ 2 =1, and |c# L | 2 
= |cl^| 2 - For a given electric field E : the amplitude of 
\ c lh\ 2 (\chl\ 2 ) shows the probability another mode con- 
tributes to the phonon wave function and is an indication 
of the mixing effect. 

Figures 1(c) and 1(d) show the calculated \cll\ 2 
and \clh\ 2 as a function of E, respectively. When 
E < 2 V/nm, the mixed component \clh\ 2 increases 
parabolically and a polynomial fitting yields | clh(E)\ 2 
= 5.6 xl(T 4 £ 2 with E in V/nm. When E > 2 V/nm, 
\clh\ 2 increases linearly with respect to E. The slope 
is found to be 2.7xl0 -3 . Overall, the mixed compo- 
nent is within 1.2% with E up to 5 V/nm, implying a 
small mixing effects induced by EEF on the in-plane op- 
tical modes. In Ref. [15], Ando and Koshino analyzed 
the phonon Green's functions of in-plane modes in BLG 
and reported that these modes are strongly mixed by 
the asymmetrical potential difference arising from bot- 
tom gate doping and external electric field [15| . To clarify 
the effects of the charge doping and the external electric 
fields, we performed further calculations to investigate 
the phonon eigenvectors for a charged BLG under an ex- 
ternal electric field. Our results show that the phonon 
mode eigenvectors are still only slightly mixed (smaller 
than 1%) by the external electric field along with the 
charge doping. 

Similar analyses can be applied to the out-of-plane 



where e m k is the energy of an electronic state \mk) with 
crystal momentum k and band index m, /(e m k) the cor- 
responding Fermi occupation, and r] is a positive infinites- 
imal. A factor of 2 accounts for electron spin. For 
a given mode uo = co>o, the phonon linewidth is 7 = 
— 2Im(II qi/ (cJo)) and the phonon frequency shift is Au 
= ^[Re(II qi ,(a;o)|£; F - n qiy (cj )U F =o]- 

The EPC matrix elements in Eq. (1) are calculated 
using the frozen-phonon approach developed in our pre- 
vious work [29|. We have considered the cases with and 
without the mixing effects. Only a negligible difference 
has been noticed since the mixed amplitude is within 
1.2% for the range of the electric fields considered. Be- 
low we only show the results with mixing effects included. 
Calculations of the phonon self-energy have been car- 
ried out on a dense 101x101 k-grid within a minizone 
(0.2x0.2) enclosing the BZ corner K{K') in the recipro- 
cal space. This is equivalent to 500x500 &-grid sampling 
in the whole Brillouin zone. By changing the Fermi level 
E F in Eq. (1), the dependence of 7 and Auj on differ- 



modes, the low-branch A\ 9 and the high-branch A2 U 
modes. In Figs. 2(c) and 2(d), the amplitude of the mixed 
components \cll\ 2 and \clh\ 2 are presented as a func- 
tion of respectively. When E < 2 V/nm, the mixed 
component \clh\ 2 increases parabolically as | clh(E)\ 2 
= 5.9xl0- 3 £ 2 with E in V/nm. When E > 2 V/nm, 
\clh\ 2 increases linearly with a slope of 2.2x10" 2 . From 
Figs. 2(c) and 2(d), EEF has dramatic effects on the 
mixing of the two out-of-plane modes; the mixing com- 
ponents reach 10% as E = 5.5 V/nm. This is nearly 
an order of magnitude larger than that for the in-plane 
modes. It was previously assumed that the electric field 
has negligible effects on the phonon eigenmodes 0. In 
contrast, as demonstrated here, the out-of-plane phonon 
modes in layered systems could be significantly altered 
by EEF. 

We then turn to the electron-phonon coupling (EPC) 
of the two in-plane modes. EPC plays a key role in under- 
standing many phenomena @, 0,M 0, Q3 , especially the 
Raman frequency shift and broadening. The phonon self- 
energy II qi/ (co>) of a phonon with wave vector q, branch 
index z/, and frequency Co> qzy provides the renormalization 
and the damping of that phonon due to the interaction 
with other elementary excitations. Following the Migdal 
approximation, the self-energy induced by the EPC in 
BLG reads 0, El|: 



(1) 

I 

ent doping levels can be investigated, assuming the EPC 
matrix elements are unchanged. This approximation is 
justified by the small dependence of the EPC matrix el- 
ements on doping for the T phonon modes in graphene 
[33| . For all the linewidths calculated below, a parameter 
of r] = 5 meV has been used. 

Figures [3^a) and [3jb) show the calculated linewidth 
7 for the low-branch (E g ) and high-branch (E u ) modes 
as a function of the Fermi level Ep and the electric field 
strength E, respectively. When the system is neutral (Ep 
=0), the linewidth 7 of the E g mode increases from 9.0 
cm -1 to 10 cm -1 as E increases from zero to 1 V/nm. 
It increases to the maximum of 30 cm -1 when E = 3 
V/nm. After E > 3 V/nm, 7 decreases to zero rapidly. 
As shown in Fig. 01(b), the largest 7 at E =3 V/nm is 
due to the fact that the electronic transition amplitude 
is the highest when the field-induced band gap ~ 0.2 eV, 
close to the phonon energy. In contrast, when Ep =0, the 
high branch E u mode exhibits nearly negligible linewidth 
for all the electric fields. This result clearly shows that 
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only the low-branch E g mode can possibly exhibit a large 
linewidth as observed in previous IR measurement [9j, 
and thus resolved the discrepancy on the mode assign- 
ment 0, [To[ . The calculated phonon linewidth without 
field reproduces the features found in the previous DFT 
[35| and analytical calculations [14] and agree reasonably 
well with the experimental data [7[. The evolution of 7 
for the Eg mode as a function of E also agrees quantita- 
tively with that by Ando and Koshino [15| . 

For a given field E, the linewidths of the two modes 
change dramatically with the doping level Ep. This can 
be understood from the selective coupling between the 
phonon modes and electronic bands. In the low dop- 
ing regime with \Ep\ < huo/2 ~ 0.1 eV, the low-branch 
mode can be in a resonant coupling with the electron- 
hole pair from the top valence band (vl) and the bottom 
conduction band (cl), as shown in Fig. 4. As a result, 
the linewidth of the low-branch mode is a constant within 
this doping range, as depicted in Fig. [3fa). Beyond this 
range, the linewidth of E g mode becomes nearly zero, 
while 7 of the E u mode increases. This is because E u 
mode couples with the two valence bands (vl and v2) 
and conduction bands (cl and c2), respectively, which 
has been discussed in detail in our previous work |3lj |. 
This feature is not captured by Ref. [15j because of their 
simplified band structure of BLG. In view of the slight 
mixing effect on the phonon mode, we conclude that it 
is mainly the band structure effect that plays a domi- 
nant role in the tunable interplay between electrons and 
phonons in BLG induced by EEF and Ep. 

The corresponding frequency shifts Auj are presented 
in Figs. [3fc) and [3fd), respectively. In particular, the 
low-branch E g mode exhibits a frequency increase (hard- 
ening) when the \Ep\ > huo/2. In contrast, the high- 
branch E u mode softens with the increase of \Ep\- As a 
result, two distinguishable peaks will appear in the Ra- 
man or IR spectra of gated BLG. Since the frequency of 
the E u mode (1592 cm -1 ) is higher than the E g mode 
(1587 cm -1 ), one expects a crossing of the two branches 
as \Ep\ increases. This is in agreement with the Raman 
data by Yan et al. [loj. From Figs. [3{ c) and^d), for a 
higher E, the E g mode can be tuned by \Ep\ and hard- 
ens even faster, implying an interesting interplay between 
EEF, electrons and phonons in BLG. 

In summary, we found that the electric field only 
slightly mixes the two in-plane modes, while inducing 
a relatively large change on the out-of-plane modes. Be- 
cause of the broken symmetry, the two in-plane modes 
become Raman and IR active, and are detectable by 
both Raman and IR spectroscopy. The delicate electron- 
phonon coupling can be tuned through doping and an ex- 
ternal gate field. It is mainly the band structure of BLG, 
tunable by external field and charge carrier densities, 
that dominates the dependence of the phonon linewidth 
and frequency shift as a function of Ep and E. Our find- 
ings have revealed the external field effects on the phonon 




FIG. 3: (Color online) Calculated phonon linewidth 7 of the 
(a) low-branch and (b) high-branch in-plane modes as a func- 
tion of the Fermi level Ef as well as the electric field E (in 
V/nm). The corresponding frequency shifts Auj are shown 
in (c) and (d), respectively. The neutrality point has been 
shifted to zero. 
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FIG. 4: (Color online) Band structure of BLG under (a) zero 
electric field and (b) electric field of E — 3 V/nm. As indi- 
cated in (a), the low-branch E g mode couples with transitions 
I (vl-cl), II (v2-c2), while the high-branch E u mode couples 
with transitions of III, IV, V, and VI. The neutrality point 
has been shifted to zero. 

frequency renormalization observed in experiment, and 
may also be useful to study the field effects on other lay- 
ered systems such as multilayer graphene. 
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